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We theoretically study the propagation through a resonant absorbing medium of a time-dependent 
perturbation modulating the amplitude of a continuous wave (cw). Modeling the medium as a two- 
level system and linearizing the Maxwell-Bloch equations for the perturbation, we establish an exact 
analytical expression of the transfer function relating the Fourier transforms of the incident and 
transmitted perturbations. It directly gives the gain and the phase shift undergone in the medium 
by a harmonic modulation. For the case of a pulse modulation, it enables us to determine the 
transmission time of the pulse center-of-mass (group delay) , evidencing the relative contributions of 
the coherent and incoherent (population) relaxation. We show that the group delay has a negative 
value (fast light) fixed by the coherent effects when the cw intensity is small compared to the 
saturation intensity and becomes positive (slow light) when this intensity increases, before attaining 
a maximum that cannot exceed the population relaxation time. The analytical results are completed 
by numerical determinations of the shape of the transmitted pulses in the different regimes. 

PACS numbers: 42.25.Bs, 42.50.Hz, 42.50 Gy 
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I. INTRODUCTION 



Convincing demonstrations of fast light in a linear ho- 
mogeneous medium were performed in the 1980s by ex- 
ploiting the steep anomalous dispersion associated with 
a well-isolated, narrow and strong absorption line of the 
medium, leading to negative values of the group veloc- 
ity [ll, Bl- Ideally bell-shaped pulses can then propagate 
moderately distorted in such a way that the maximum 
of the transmitted pulse occurs before that of the inci- 
dent pulse. Time-advances exceeding 0.5 times the half- 
duration at half-maximum of the pulse-envelope have 
been so evidenced in experiments involving a true shape 
detection of the latter Q- See also [J, IJ. More generally, 
due to the causality principle (originating the Kramers- 
Kronig relations), fast light is expected every time that 
the carrier frequency of the pulses coincides with the bot- 
tom of a well-marked dip in the medium transmission 
and the pulse distortion resulting from the first order 
variations of the transmission and of the group velocity 
versus frequency cancels when this coincidence is exact 
[5|. The dip in the transmission may be natural or cre- 
ated in various arrangements involving electromagneti- 
cally induced absorption ^, 7], stimulated Raman ^8| or 
Brillouin Q scattering, etc. The literature on fast light 
is abundant (for reviews, see, e.g., J10l-[l3]) but it seems 
that the fractional pulse-advance with moderate distor- 
tion reported in [3| has not been overtaken. This is eas- 
ily explained by remarking that, since the transmission 
dip at the carrier frequency must be well pronounced, 
there are spectral regions where the overall transmission 
of the system is much larger. For obvious reasons of 
noise (no matter its origin) , instability and hypersensitiv- 
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ity to small defects in the incident pulse-shape, the usable 
transmission-dynamics cannot be too large and this lim- 
its the observable fast-light effects in any linear system 
[l4| . The dynamics actually involved in the fast-light 
experiments does not exceed 45 dB and is often much 
lower. The situation is quite different in the slow-light 
experiments where the medium transmission is maximal 
at the carrier frequency and can be very low on both sides 
of this frequency, with transmission-dynamics exceeding 
1000 dB 0,[i||. 

In the present article, we study the propagation in an 
absorbing medium of a continuous wave (cw) whose am- 
plitude is pulse-modulated with a low m,odulation index. 
The medium is modeled as a homogeneously-broadened 
two-level system (TLS) [13| with a resonance frequency 
loq and a relaxation time Ti (T2) for the population dif- 
ference (the coherence). We assume that the cw is on 
exact resonance and we are interested in the propaga- 
tion of the pulsed part of the wave (the pulse) whose 
optical spectrum is centered on ujq and maximum at this 
frequency. The pulse may thus be seen as a continu- 
ous superposition of fields symmetrically detuned from 
the cw field (sidebands) that probe the changes of the 
TLS properties induced by the cw. The propagation in a 
strongly driven TLS of two symmetrically detuned fields 
is a basic problem in nonlinear optics [1^ 13 • When the 
sidebands are associated with an amplitude modulation, 
it has been experimentally evidenced [20| and theoreti- 
cally demonstrated [isl . [2l| that the transmission of the 
modulation presents a maximum on resonance for strong 
enough cw intensities. According to our general anal- 
ysis, the pulse propagation may be slow in such condi- 
tions. On the other hand, when the cw intensity is small 
compared to the saturation intensity, the TLS behaves 
linearly both for the cw and the pulse, the propagation 
of which is not modified by the presence of the cw (fast 
light regime) . Our main purpose is to analyze how takes 



place the transition from a fast to a slow light regime 
when the cw intensity increases. Our calculations take 
full account of the depletion of the cw intensity during 
the propagation and are made for arbitrary values of the 
relaxation times Ti and T2. The corresponding results 
(mainly analytical) are new and it seems that the pos- 
sibility of switching from fast to slow light in a TLS by 
simply adding a resonant cw to the pulses has not been 
previously considered. 

The arrangement of our paper is as follows. In Sec. |lTl 
we specify the system under consideration and we estab- 
lish an exact analytical expression of its transfer function 
for the time-dependent part of the field envelope. This 
expression directly gives the gain and the phase shift un- 
dergone in the medium by a harmonic modulation, ex- 
amined in Sec lIIII The propagation of pulses is studied in 
Sec. IIVI with a special attention paid to the transmission 
delay of their center-of-mass (identified to the group de- 
lay). The analytical results are completed by numerical 
determination of the envelope of the transmitted pulses. 
We finally conclude in Sec. |V]by summarizing the main 
results. 



II. TRANSFER FUNCTION OF THE MEDIUM 

As above mentioned we consider the case where the 
cw is exactly resonant. The amplitude modulation is 
then an eigenmode of modulation in the sense that that 
it is conserved during the propagation [23 - t24| . We de- 
note £ the thickness of the medium, a its unsaturated 
absorption coefficient on resonance for the intensity and 
z (0 < z < i) the direction of propagation of the wave 
assumed to be plane and polarized along the x-axis. We 
write the x-component of the electric field as : 



E^{z,t) = Re e"^<'*E{z,t) 



(1) 



where E{z,t) is the slowly varying field-envelope. As 
in all the following, t is a local time, that is the real 
time minus ngz/c, c being the velocity of light in vacuum 
and no the refractive-index around loq of the eventual 
host medium. Denoting /i the dipole matrix element for 



the transition (chosen real), R{z,t) 



tJ,E{z,t) 



the Rabi 



frequency, N{z,t) the population difference per volume 
unit (A^o its value at equilibrium) and P{z,t) the enve- 
lope of the electric polarization induced in the medium, 
it is convenient to introduce the dimensionless quantities 

Note that all these quantities are real for the resonant 
case considered here and that / = _E^ is the intensity 
normalized to the saturation intensity. In such condi- 
tions, the Maxwell-Bloch (MB) equations governing the 
system evolution take the simple form : 



dE 
dz 



-P 



(2) 



dP 
T2^^DE~P 
ot 



dD 

Ti— = -PE + {1 - D) 



(3) 



(4) 



These equations are easily solved when only the cw is 
present. In the following we denote the corresponding 
values of E, P, D and / by the index cw. Combining 
Eqs. © and (g]), we find Dcwiz) = 1/[1 + Icwiz)] and 
Pcwiz) ~ Ecn]{z)/ [I + Icw{z)]. Replacing this result in 
Eq. dll), we get : 



d 

— [Icw{z) + ln/c^,(z) + az]=0 



In 



InL, 



In I,, 



ai 



(5) 



(6) 



where lout and Im are short hand notations of Icwi^) 
and /ct«(0). When a small time-dependent perturbation 
SEin (t) (also real for an amplitude modulation) is added 
to the constant amplitude E^ = y/Tin of the field of the 
incident cw, we search for solutions of the MB equations 
under the form E{z,t) = Ecw{z) + 5E{z,t), D{z,t) = 
Dcyjiz) + 5D{z,t) andP(z,t) = Pcyjiz)+SP{z,t). Taking 
into account the cw solution, we get at the first order of 
perturbation 



diSE) 



a 



5P 



T2^^ = E,u:SD + D,,,SE - 5P 
at 



Ti^^ = -E^^SP - Pc^^SE - SD 



dt 



(7) 



(8) 



(9) 



These linear equations can be solved by introducing the 
Fourier transforms /S.E{z, O), AP(z, Vt) and AD(z, Vt) of 
6E{z,t), SP{z,t) and 6D{z,t) Hg. From Eqs §1 and 
®, we get : 



AP{z,n) 



[1 + iilTi 



IcUz)]AE{z,n) 



[1 + Ic^iz)] [(1 + inTi) (1 + 1^X2) + Icw{z)] 

(10) 
Replacing this result in Eq.(I7]) and integrating on z, 
we finally obtain the transfer function H{Vt) relat- 
ing A£'otjt(r2) and A£'i„(r2), the Fourier transforms of 
5Eout{t) = SE{i,t) and 6Ein{t). It reads as : 



H{n) = exp j I f (z, n) dz 



(11) 



where 



f{z,n) 



a 



Icn,{z) - [l + iVlTi) 



2 [1 + /e^(z)] [(1 + im^) (1 + inT2) + Icn,{z)] 

(12) 



is the complex gain factor at the abscissa z. By means 
of Eq. ([5]) , the integration on z in Eq. pT|) can be trans- 
formed in an integration on lew ■ We get : 



ff (17) = exp 



and finally 



H{?L) = exp 



1 + iVLTx ~ Ic 



hw [(1 + i^Ti) (1 + inT2) + lew] 



dir 



(13) 



l + «^, ( hn + {l + i^Ti){l + iVlT2) 



1 + iVLT2 



■In 



h(l + zfiTi)(l 
In yj hn/ lout 

{i + inT2) 



inT2) 
(14) 



where /oMt and /„ are related by Eq.®. This expres- 
sion of H{n) is the central result of our paper. In the 
limit lin <C 1, lout ~ hn exp(— a£) and H{fl) is reduced 
to exp(— „,. "f,^ , ) which, as expected, is the transfer 
function of a linear medium with a Lorentzian absorp- 
tion line [5|. Quite generally, we see that the poles 
of if (ri) are in the half-plane Im (fi) > and that 
H (— SI) = H* (n). This ensures that the impulse re- 
sponse h{t), inverse Fourier transform of H{n,), cancels 
for i < and is real, confirming that the system is causal 
p5] and that the amplitude modulation is an eigenmode 
of modulation. At the first order of perturbation consid- 
ered here, we finally remark that 



I{z,t) = [Ecwiz) + 5Eiz,t)f « Icwiz)+26E{z,t)^Icw{z) 

(15) 
This means in particular that a quadratic detection, as 
currently used in optics, will deliver a time dependent 
signal SIout{t) proportional to 5Eout{t) and that the in- 
stantaneous modulation indices for the intensity and for 



the field are such that 



SI(z,t) _ r,SE(z,t) ■ 



»w 



= 2 



Ecru{z) 



in every point. 



III. TRANSMISSION OF A HARMONIC 
MODULATION 
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of the output amplitude-modulation with respect to the 
input one. The gain coefficient g{z, fl) = Re [f{z, SI)] and 
the phase-shift coefficient ip{z, il) = Im [f{z, i7)] are the 
corresponding quantities per unit length. Due to the de- 
pletion of the cw intensity lew during the propagation, 
these latter quantities depend on z. Note that the expres- 
sion oi g{z, H.) derived from Ea. ([T^ is consistent with the 
results given in [1^ [21| , also as in [20] for the particular 
case Ti = T2. Figures [2] and [3] show a set of profiles g{fl) 
and (/3(r2) obtained for T2 = 2Ti (purely radiative relax- 
ation) with various cw intensities. For lew ^ 1 (fully 
linear case), we obviously retrieve the gain and disper- 
sion profiles g(f]) = -^^^^^ and ^{n) = ^^^^^ , 

associated with a Lorentzian absorption line. 




Frequency 



Figure 2: Gain coefficient g(0.) in a/2 units as a function of 
Q in l/Ti units for T2 = 2T\. The different curves are labeled 
by the corresponding value of the normalized intensity lew ■ 
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Figure 3: Same as Fig[2] for the phase-shift coefficient tf(Q,) 
in a/2 units. 



Figure 1: The driving wave (cw) at ujq can modify the prop- 
agation of the pair of sidebands at ujo ±0.. 

When the amplitude-modulation is harmonic of fre- 
quency fi (f2 > 0) , the optical field has two sidebands at 
wq ± i7 that both act as probes of the system driven at wq 
(FiglH). The modulus G{Vt) and the argument $(f2) of 
H{il) are then respectively the gain and the phase-shift 



The gain coefficient (Figl5]) keeps negative as long as 
lew < 1. It cancels at fi = when lew = 1 (absolute 
cw intensity equal to the saturation intensity). Beyond, 
it becomes positive in the frequency range Q, < Qg = 

Vi^cw ~ 1) / {IcwTiT2 + T^), attaining an absolute max- 
imum a/16 for lew — 3. When lew ^ 1, Slg « Rew 
, the Rabi frequency associated with the cw field. We 
also note that, g{fl) being an even function, g{0) = 



"/, """.^X is always an extremum of the gain-coefficient. 

This extremum is flat for lew = b ± ^/W^^T^/Ti with 
b ^ 3/2 + T1/T2 + T2/2T1, that is lew = 3±V7 when 
T2 = 2Ti . Despite some common points, the profiles g(ri) 
significantly differ from those obtained by using a probe 
independent of the driving field [isl . [26l - [29| . They are a 
bit simpler, probably because the amplitude modulation 
is an eigenmode of modulation (contrary to the single 
sideband modulation). On the other hand, the refractive 
index experienced by the optical probe field, well defined 
for the single sideband case [30|, is not defined for the 
case of an amplitude modulation where the probe field 
has two frequency-components ed loq ± D, that are not 
independent. 

The phase-shift coefficient (p{ft) (Fig. [3]) has a max- 
imum (fmax — ol/A (phase advance) ai Q. — \/T2 for 
lew ^ 1 and an absolute minimum Lpmin ~ —0.078 a 
(phase lag) at J7 « 1.14/Ti for lew ~ 3.3. On the other 



hand, the slope 



in 



n=o 



at the origin is also maximal for 



lew "^ 1- This slope is half- maximum, cancels and at- 
tains its (negative) minimum, for lew respectively equal 
to 0.135, 1/2 and 5/4. 
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Figure 4: Gain G{Q) as a function of f2 in l/Ti units for 
T2 = 2Ti and a£ — 10. The different curves are labeled by 
the corresponding value of the normalized intensity Ii„ . 



Representative profiles of the (overall) gain G{il) and 
phase-shift $($7) are shown Figl4]and FigOfor a£ — 10, 
with various intensities /„ chosen according to crite- 
ria analogue to those used Fig 15] and Fig 121 The op- 
tical thickness retained is of the order of that actually 

used in [3|- When /, 
and $(17) = 



< 1, Gin) = exp 



_ ainT-i 



Except for this fully linear 



2(l+n2T|)' 

case, Eqs. dlT]) and ^2^ show that G{D,) and $(12) are 
strongly affected by the depletion of the cw intensity in 
the medium, often not taken into account in the litera- 
ture. This point is illustrated FiglS] where we give the 
profile 0(17) = exp [ig{il.)] which would be obtained by 
neglecting the cw depletion, i.e. by taking Iew{z) — /,„ 
everywhere. The depletion obviously forces one to use 
larger incident cw intensities to obtain similar gain pro- 
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Figure 5: Same as Fig|4]for the phase-shift $($!). The maxi- 
mum and minimum phase-shift are respectively af /4 = 2.5 Rd 
at n = IjT'i for /i„ < 1 and -0.56 Rd at fi = 1.68/Ti for 
/,„ ^ 11.5. 



files and the maximum gain is smaller. The difference 
is more important when the optical thickness at is very 
large. From Eqs. Q and p4p . we find that the maximum 
gain then tends to \/al/2 whereas it equals exp [a£/16] 
when the depletion is neglected. There are obviously sim- 
ilar effects on the phase-shift $(51). 
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Figure 6: Same as Fig|4] when the depletion of the driving 
wave is not taken into account. 

The previous results take a very simple form in the 
limit case where T2 is negligible compared to Ti . The MB 
equations are then reduced to the rate equations |17| that 
provide a good description of th e propa gation of light in 
numbers of saturable absorbers |3ll - l36| . Putting r2 ~ 
in Eqs. ^ and dH]), we get 



f{z.n) 



ale 



H{n) 



2{1+Iew) [l+I^^f \l+iJ^^ 

+ lin + i^Ti 



^out 
^ in. 



1 



iQTi 



Apart from the term 



- 2(1+1 — ) independent of fl, 
the complex gain factor f{z,fl) is identical to that 



of a Lorentzian gain line of half-width at half maxi- 
mum (1 -l- lew) /Ti with a gain-coefficient on resonance 
dlcw/ (1 + lew) for the amplitude. During the propaga- 
tion in the medium, the modulation index is magnified 
whereas the phase of the modulation is delayed [35. - .38| . 
The magnification K{U) of the modulation index and the 
phase lag i(J7) — — <J>(J7) for both the amplitude and the 
intensity modulations are easily deduced from Eq. (|17|) . 
We get : 



K{n) 



G{n)^l-^ 

^out 



' {i+hn)'+n^n 



L[il) = tan 



r!Ti (/,, 



[1 + Ua) {I + lout) + n^Tl 



(18) 



(19) 



The phase lag, always positive for this "inco- 
herent" case, attains its maximum Lm.n.T. = 
-1 



tan 



{hn - lout) l2yJ{l+hn){l+Iout) 



for 



VlTi = ^(H-/j„) (1 + lout)- We remark that 
Lmax < 7r/2, the upper limit being approached when 
lin ^ 1 and lout '^ 1) that is for a€ — >■ oo. Consequently 
the time delay of the output modulation can never 
exceed a quarter of the modulation period T — 27r/J7 
(about an eighth for ai = 10). Since the work reported 
in [33, numerous slow-light experiments performed in 
saturable media have been analyzed by invoking hole 
burning via coherent population oscillations (CPO), 
resulting in a reduction of the group velocity. In fact, a 
more direct analysis is provided by the basic model of 
saturable absorber [33J and, as shown in detail in [36j, 
the signals observed in most CPO experiments can be 
perfectly reproduced by means of Eqs. (fT8|) and ([T9| . 
eventually extended to take into account inhomogeneous 
effects. 



emphasize that, due to the cw depletion, Vg is not uni- 
form. From Eq. (|14p , we get Tg — ti + T2 with 



ri=Ti 



1 



1 



T2 = T2 



lout + 1 lin + 1 



-In 




(20) 



(21) 



Ti (T2), proportional to Ti (T2), may be considered as 
the contribution to the group delay of the incoherent 
(coherent) effects. We see that the incoherent or popula- 
tion effects always lead to a delay whereas the coherent 
ones mainly lead to a much larger time-advance (FiglT]), 
at least when T2 is comparable to Ti (the coherent ef- 
fects obviously disappear in a saturable absorber where 
T2 « 0). 
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Figure 7: Normalized coherent and incoherent contributions 
to the group delay as a function of the incident driving in- 
tensity lin, depending whether the cw depletion (q£ = 10) is 
taken into account or not. 



IV. PULSE PROPAGATION 

Strictly speaking a harmonic modulation does not con- 
tain any information and, e.g., a time-delay of T/8 as con- 
sidered at the end of the previous section can also be seen 
as a time-advance of IT 1%. Unambiguous demonstrations 
of fast or slow light require to use pulses of finite duration 
and energy. An important parameter to characterize the 
propagation is then the transmission-time of the center- 
of-mass of the pulse envelope. Following the use in sig- 
nal theory [25|, we define the center-of-mass of a signal 
y(t) of normalized area as /_ ty{t)dt. A direct appli- 
cation of the moment theorem [23 shows then that the 
transmission-time of the center-of-mass of the pulse en- 
velope is equal to the group delay Tg = —^ \n=Q , what- 
ever the pulse distortion may be |5j]. On the other hand 
— ^ |n=o appears as the transmission delay or group de- 
lay per unit length, equal in our local time picture to 
1/vg — no/c ~ 1/vg, where Vg is the group velocity. We 



Before commenting further on Eqs. ([20| and (|2T|). con- 
sider the simpler results obtained when the depletion of 
the driving wave is neglected {lew — lin everywhere). 
From Eq.(|12p. we then get : 



n^Ti 



alL, 



T2 



'T: 



(1 + hn)" 



ai{l-hn) 
'2(l + /„)^ 



(22) 



(23) 



Ti/Ti, negligible for /i„ ^ 1, is maximum for /.;„ ~ 1/2 
whereas T2/T2 starts from the large negative value —a£/2 
before attaining a small positive maximum for I in = 2. 
The total group delay cancels for lin = T2/ (T2 + 2Ti) 
and is maximum for /i„ = (2r2 + Ti) / {T2 + 2Ti) that is 
respectively for /„ — 1/2 and lin = 5/4 when the relax- 
ation is purely radiative (see Figl3]). As expected, taking 
into account the depletion of the cw intensity does not 
modify the delays for lin ^ 1 but shifts the curves Ti{Iin) 



and T2{Iin) to larger /i„, dramatically spreads them and 
reduces the amplitude of their maximum (Fig|7]). Equa- 
tion ([20|) shows that ri cannot overtake Ti, this value 
being approached when /„ 3> 1 and lout ^ 1, that is for 
ai extremely large J35| . An estimate of the maximum of 
Ti for large but realistic values of a(. can be obtain by 
anticipating that it is attained for lout ^ 1- From Eq.(p|) 
we get then lout ~ hn exp {hn - at) and 



Tl«Tl 



1 



lin exp [lin - ai) + 1 lin + 1 



(24) 



The maximum is obtained for (/j„ + 1) ~ 

exp (a£ -/„,). For al = 10 (FigO, we find /»„ w 4.8 
and, putting this value in Ea. (l24)) . max(ri/ri) « 0.80 
in perfect agreement with the exact value. By the same 
method we find max(Ti/ri) « 0.995 for at = 200. An 
upper bound to the total group delay can be obtained by 
using the method of the Lagrange multipliers, perfectly 
adapted to the search of a maximum under constraint, 
here fixed by Eq.®. Following the procedure exactly as 
exposed in the original work of Lagrange [40|, we find 
an extra relation between /,„ and lout at the maximum 



of Tg , namely 



^out 



1 



2Ti + Ta (1 - l/Iout) 

2T^+T2{l-l/hn) 



(25) 



A numerical exploration shows that the upper bound to 
the group delay is attained when /i„ is extremely large 
whereas lout keeps finite. Equation (I25|) takes then the 
form 

(1 + loutf ~ I?n [lout - T2/ (2ri + T2)] (26) 

which can be verified if and only if lout ^ T2/ (2Ti +T2). 
Injecting this result in the general expression of Tg, we 
finally find : 



T2 
sup [max (Tg)] = Ti + — - r2 In 




(27) 



When T2 « (saturable absorber) , we obviously retrieve 
the upper bound Ti for Tg whereas this upper bound 
is (2 — ln3)ri « G.9OT1 when the relaxation is purely 
radiative (T2 = 2Ti). For the latter case max(Tg/Ti) is 
only 0.68 for a£ = 10, significantly below its upper limit, 
and raises to 0.89 for aS = 200. 

As before mentioned, the group delay Tg is the 
transmission-delay of the pulse center-of-mass, as large 
as the pulse distortion may be. On the other hand, 
the envelope SEoutit) of the transmitted pulse is sim- 
ply the inverse Fourier transform of _ff(51)A£'i„(f7). 
When the duration Tp of the incident pulse is long 
enough, A£'j„(f2) is concentrated around 17 = where 
H {n) Ki y/ Iin/Ioutexp{~inTg). We then get 5Eout{t) = 
-\/Iin/IoutSEin {t — Tg). The whole envelope is multiplied 

by \/ lin I lout and time shifted by r^, without any distor- 
tion. Strictly speaking, this solution is only valid for a 



very long pulse and is not really interesting insofar as the 
time delay is then negligible compared to the pulse dura- 
tion. For finite Tp, there is always some pulse-distortion 
whose importance depends in particular on the trans- 
mission dynamics (as defined in the introduction). If, as 
usual, the incident pulse is bell-shaped and symmetric 
with a maximum at i = 0, the envelope of the trans- 
mitted pulse will be dissymmetric with a maximum at 
a time r^ such that Tm/Tg < 1, its center-of-mass keep- 
ing exactly at the time Tg no matter the pulse duration. 
The challenge in the slow or fast light experiments is to 
obtain a fractional delay or advance |t„i| /Tp as large as 
possible, with moderate distortion, t^ and Tg are then 
comparable. The following figures are obtained by us- 
ing standard techniques of Fast Fourier Transform. The 
incident pulse is Gaussian and Tp is its half-duration at 
half-maximum. 

FigureslSl [9land[T0lshow typical shapes of the transmit- 
ted pulse for T2 — 2Ti (radiative relaxation) and a£ = 10. 
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Figure 8: Normalized pulse envelopes obtained for a£ = 10 
and T2 — 2T\, when the group advance —Tg is maximum 
(solid line) and half-maximum (dashed line), r^ is chosen such 
that the fractional advance of the pulse maximum —Tm/Tp is 
significant whereas the distortion remains moderate {rp/Ti = 
12.5). For the first case, Tm/Tg — 0.88 and —Tm/Tp = 0.71. 
For the second case TmJTg = 0.91 and —Tm/Tp = 0.37. The 
envelope of the incident pulse is given for reference (dotted 
line) . 



Figures [5] and [S] are obtained for lin — 0.001 and 
lin — 1-9 leading respectively to a group advance —Tg 
nearly equals to its maximum a£T2/2 and the half of 
this value. The corresponding transmission dynamics is 
respectively 43 dB and 26 dB. On Fig[51 Tp has been 
chosen in order that the fractional advance —Tm/Tp is 
significant (respectively 0.71 and 0.37) and the distor- 
tion keeps moderate (jm/Tg respectively equals 0.88 and 
0.91). When the pulse duration is shortened, the frac- 
tional advance increases but not as much as one could 
expect (respectively up to 1.1 and 0.68) and this is paid 
by a dramatic distortion of the transmitted pulse (Figl9]) . 

Figure [10] is obtained for /„ sa 8.0, leading to the 
largest normalized group delay (Tg/T\ k, 0.68). The cor- 
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Figure 9: Same as Fig. 8 when the pulse durations are those 
maximizing the fractional advance —Tm/fp . For the solid 
line Tp/Ti = 4.1 , Tm/rg = 0.44 and -Tm/rp = 1.1. For the 
dashed line Tp/Ti = 2.8 , Tm/Tg = 0.38 and —Tm/rp — 0.68. 
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Figure 11: Same as Fig [10] when T2 ~ Q (saturable absorber) 
For the solid line Tp/Ti = 2.0 , Tm/Tg = 0.79 and TmJTp 
0.32. For the dashed line Tp/Ti = 0.44 

TmlTp = 0.57. 



/Tg — 0.31 and 



CD 
Q. 
O 
CD 
> 
C 
LU 
(D 

jn 
a. 



1.0- 


/ 


/■■ 

// 

.■■■ / 
// 
I 




\ \ 


0.5- 


/ 








0.0- 


1 1 




1 


1 1 1' 



Time (in Xp-units) 



Figure 10: Normalized pulse envelopes obtained for a£ = 10 
and T2 = 2T\, when the group delay Tg is maximum. For the 
solid line Tp/Ti = 2.27 , Tm/Tg = 0.77 and t^/tp = 0.23. For 
the dashed line Tp/T\ = 1.44 , Tm/Tg = 0.56 and Tm/Tp = 
0.27. The envelope of the incident pulse is given for reference 
(dotted line). 



responding transmission dynamics is small (7dB). This 
explains in part that, even when the pulse duration is 
optimized (dashed line), the fractional delay T^n/Tp does 
not exceed 0.27. Another reason is that the coherent re- 
laxation negatively contributes to the group delay (see 
FigED. 



For a saturable absorber (T2 « 0) of same optical 
thickness «£, the maximal group delay, now attained for 
lin ~ 4.8, is upgraded to O.8OT1 (see FigjT]) whereas the 
transmission dynamics attains 15 dB. Figure [Til shows 
the pulse envelopes obtained for two different pulse du- 
rations. As expected, the maximum fractional delay, ob- 
tained for Tp K, 0.44ri (dashed line), is significantly larger 
than for the case of radiative relaxation (more than two 
times larger) . We however remark that the corresponding 
fall of the pulse is considerably lengthened. An extensive 
study of the pulse distortion in this particular system can 
be found in [35]. 



V. CONCLUSION 

In their letter appeared in 1982 and soberly entitled 
"Linear Pulse Propagation in an Absorbing Medium", 
Chu and Wong [1| paved the way for the now-called fast- 
light experiments. We have examined in the present ar- 
ticle how the propagation is modified when the absorb- 
ing medium is driven by a cw. Modeling the absorbing 
medium as an ensemble of two-level atoms, we have more 
specifically considered the case where both pulses and 
cw are on exact resonance, a condition simply realized 
by pulse-modulating the cw amplitude, with a low mod- 
ulation index. This particular arrangement eliminates 
the pulse-distortion associated with the first order vari- 
ations of the absorption and of the group delay versus 
frequency. The basic result of our paper is the exact an- 
alytical expression of the transfer function relating the 
Fourier transforms of the incident and transmitted mod- 
ulations for arbitrary values of the coherent and incoher- 
ent relaxation times (Ea ll4p . It shows the importance of 
the effects resulting from the depletion of the cw inten- 
sity along the propagation. When the modulation is only 
harmonic (as in numerous experiments), it directly gives 
the gain and the phase shift undergone by the modula- 
tion. They significantly depart from those obtained with 



a unique probe field (single sideband modulation). When 
the modulation is actually pulsed, the transmission de- 
lay of the pulse center-of-mass, identified to the group 
delay, is deduced from the transfer function by a simple 
calculation of derivative. A remarkable point is that the 
group delay is the sum of two terms, respectively propor- 
tional to the coherent and incoherent relaxation times. 
These two terms being mainly of opposite sign and de- 
pending differently on the cw intensity, this explains why 
the transmission delay, strongly negative when the cw in- 
tensity is low, may become (slightly) positive when the 
latter increases. Finally the numerical determination of 



the pulse shapes confirms a general property of the fast 
and slow light systems, namely that significant advances 
or delays with moderate distortion can only be obtained 
in media with a large transmission dynamics. Though 
our study is only theoretical, it is illustrated for a real- 
istic value of the optical thickness, comparable to that 
actually used in the microwave experiment reported in 
[2] . In the optical domain, suitable optical thickness and 
time-scale could probably be obtained by using an en- 
semble of cold atoms. We expect that our theoretical 
work will stimulate such an experiment. 
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